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In the present paper we present results of calculations obtained with the use of the 
theoretical method described in our preceding paper [jj and perform detail analysis 
of a-helix<->random coil transition in alanine polypeptides of different length. We 
have calculated the potential energy surfaces of polypeptides with respect to their 
twisting degrees of freedom and construct a parameter-free partition function of the 
polypeptide using the suggested method [j| . From the build up partition function we 
derive various thermodynamical characteristics for alanine polypeptides of different 
length as a function of temperature. Thus, we analyze the temperature dependence 
of the heat capacity, latent heat and helicity for alanine polypeptides consisting of 21, 
30, 40, 50 and 100 amino acids. Alternatively, we have obtained same thermodynam- 
ical characteristics from the use of molecular dynamics simulations and compared 
them with the results of the new statistical mechanics approach. The comparison 
proves the validity of the statistical mechanic approach and establishes its accuracy. 

I. INTRODUCTION 

In our preceding paper we introduced a novel and general theoretical method for 
the description of phase transitions in finite complex molecular systems. In particular, we 
have demonstrated that for polypeptide chains, i.e. chains of amino acids, one can identify 
specific twisting degrees of freedom that are responsible for the folding dynamics of these 
amino acid chains. In other words, these degrees of freedom characterize the transition from 
a chain in a random coil state, to one in an a-helix structure and vice versa. 

The essential domains of the potential energy surface (PES) of polypeptides with respect 
to these twisting degrees of freedom have been calculated and thorougly analyzed on the 
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basis of both classical molecular dynamics (MD) simulations, and ab initio methods such as 

nnn 

density functional theory (DFT) and the Hartree-Fock approach. In Refs. [1|, [2|, [3J , it was 
shown that with the PES, one can construct a partition function of a polypeptide chain from 
which it is then possible to extract all essential thermodynamical variables and properties, 
such as the heat capacity, phase transition temperature, free energy, etc. 

In this paper, we explore this further using a formalism we introduced previously B 
and apply it to a detailed analysis of the a-helix^random coil phase transition in alanine 
polypeptides of different lengths. We have chosen this system because it has been widely 



investigated both theoretically 



experimentally [21 
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during the last five decades (for review see, e.g. [ 
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and thus is perfect system for testing a novel theoretical approach. 

The theoretical studies of the helix-coil transition in polypeptides have been performed 
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both with the use of statistical mechanics methods 
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and of MD 
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Previous attempts to describe the helix-coil transition 



in polypeptide chains using the principles of of statistical mechanics were based on the 



models suggested in sixties 



These models were based on the construction of 
the polypeptide par tition function depending; on several parameters and were widely used 
in Refs. 116. 



13, 1181, 
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25l . |26| . |27[ | for the description of the helix-coil transition in 



polypeptides. 



For a comprehensive overview of the relevant work we refer to recent reviews 25|, |27J, |2_8] 



and the book 



26| . as well as to our preceding paper 



Experimentally, extensive studies o 
conducted 
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the helix-coil transition in polypeptides have been 
24(. In Ref. [2l|, the enthalpy change of an a-helix to random coil 
transition for the Ac-Y(AEAAKA) 8 F-NH 2 peptide in water was determined calorimetrically. 
The dependence of the heat capacity of the pol ype ptide on temperature was measured using 
differential scanning calorimetry. In Refs. 22l, |23|], UV resonance Raman spectroscopy was 
performed on the MABA-[A]5-[AAARA]3-ANH2 peptide. Using circular dichroism methods, 
the dependence of helicity on temperature was measured. While in Ref. 241 ]. the kinetics of 
the helix-coil transition of the 21-residue alanine polypeptide was investigated by means of 
infrared spectroscopy. 

In this work, we have calculated the PES of polyalanines of different lengths with respect 
to their twisting degrees of freedom. This was done within the framework of classical molec- 
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ular mechanics. However, to scrutinize the accuracy of these calculations, we compared 
the resultant molecular mechanics potential energy landscapes with those obtained using 
ab initio density functional theory (DFT). The comparison was only performed for alanine 
tripeptide and hexapeptide, since for larger polypeptides, the DFT calculation becomes 
increasingly computationally demanding. Hence for these larger systems, only molecular 
mechanics simulations have been used in this present work. 

The calculated PES was then used to construct a parameter-free partition function of the 
polypeptide using the statistical method we had outlined in our preceding paper This 
partition function was then used to derive various thermodynamical characteristics of alanine 
polypeptides as a function of temperature and polypeptide length. We have calculated and 
analyzed the temperature dependence of the heat capacity, latent heat and helicity of alanine 
polypeptides consisting of 21, 30, 40, 50 and 100 amino acids. We have also established a 
correspondence between our ab initio method with the results of the semiempirical approach 
of Zimm and Bragg [jj]. Thus, on the basis of our approach, we have determined the key 
parameters of the Zimm-Bragg theory that itself utilizes principles of statistical mechanics. 

Finally, we have calculated the heat capacity, latent heat and helicity of alanine polypep- 
tides using molecular dynamics and have compared the obtained results with those using our 
statistical approach. Comparison between the two methods allows us to establish the accu- 
racy of our statistical method for relatively small molecular systems, and lets us gauge the 
feasibility of extending the description to larger molecular objects for which it is especially 
essential in those cases where MD simulations are hardly possible due to computational 
limitations. 

Our paper is organized as follows. In section [III we present the final expressions obtained 
within the formalism described in our preceding paper [l| and introduce basic equations 
and the set of parameters which have been used in MD calculations. In section III 1 1 we 
present and discuss the results of computer simulations obtained with the use of developed 
theoretical method and compare then with results of MD simulations. In section IIV1 we 
draw a conclusion to this paper. 
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II. THEORETICAL METHODS 
A. Statistical model for the a-helix^random coil phase transition 

Our calculations have been performed using the statistical formalism we described pre- 
viously [lj. Here, we will only outline the basic ideas of this method and present the final 
expressions that were used in our investigation. 

Let us consider a polypeptide, consisting of n amino acids. The polypeptide can be found 
in one of its numerous isomeric states with different energies. A group of isomeric states with 
similar characteristic physical properties is called a phase state of the polypeptide. Thus, a 
regular bounded a-helix state corresponds to one phase state of the polypeptide, while all 
possible unbounded random conformations can be denoted as the random coil phase state. 

The phase transition is then a transformation of the polypeptide from one phase state to 
another, i.e. the transition from a regular a-helix to a random coil conformation. 

All thermodynamical properties of a molecular system are described by its partition 
function. The partition function of a polypeptide can be expressed as l|): 



\3N-3-if 



n— 4 



i=l 



Z = A ■ B(kT) ■ (JfcT) 

(n-3)/2 n-i-3 



(i) 



(k - l)\(n — k — 3)\ 



i\(i-l)\( k - i)!(n -k-i- 3)! 

Here the first and the third terms in the square brackets describe, respectively, the partition 
function of the polypeptide in the a-helix and the random coil phases. The second term in 
the square brackets accounts for the situation of phase co-existence. The summation in this 
term is performed up to n — 4 as the shortest a-helix has only 4 amino acids. The final 
term in the square brackets accounts for the polypeptide conformations in which a number 
of amino acids in the a-helix conformation are separated by amino acids in the random coil 
conformation. The first summation in this term goes over the separated helical fragments 
of the polypeptide, while the second summation goes over individual amino acids in the 
corresponding fragment. Polypeptide conformations with two or more helical fragments are 
energetically unfavorable. This fact will be discussed in detail further on in this paper. 
Therefore, the fourth term in the square brackets Eq. ([!]) can be omitted in the construction 
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of the partition function. 

A in Eq. (JT]) is a factor that is determined by the specific volume, momenta of inertia 
and frequencies of normal vibration modes of the polypeptide in different conformations [lj] , 
l s is the total number of the "soft" degrees of freedom in the system. B(kT) is a function 
defined in our preceding paper [lj], which describes the rotation of the side radicals in the 
polypeptide. and Z u are the contributions to the partition function from a single amino 
acid being in the bounded or unbounded states respectively. They can be written as: 



Z„ = ££e X p(-«)d^ (3) 
rj = QT £ exp p< % ,yW) ) ^ . (4) 

where k and T are the Boltzmann constant and the temperature respectively, while iV is 
the total number of atoms in the system. e^((p,ip) and e^ u \cp,ip) in Eqs. (J2J)- (J3I) are the 
potential energies of a single amino acid in the bounded and unbounded conformations 
calculated respectively versus the twisting degrees of freedom if and ip. These degrees of 
freedom are defined for each amino acid of the polypeptide except for the boundary ones 
and are described by two dihedral angels ipi and ipi (see Fig. [T|) 

Both angles are defined by four neighboring atoms in the polypeptide chain. The angle ipi 
is defined as the dihedral angle between the planes formed by the atoms (C^ — iVj — Cf ) and 
(Ni — Cf — Cj). While the angle ipi is defined as the dihedral angle between the (iVj — Cf — C'A 
and (Cf — C- — Ni + i) planes. The atoms are numbered from the NH2- terminal of the 
polypeptide and <fi and ipi take all possible values within the interval [— 180°;180°]. For 
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33j, <fi and ipi are counted 



an unambiguous definition most commonly used 
clockwise if one looks on the molecule from its NH2- terminal (see Fig{T]). 

By substituting Eqs. ([2]) and ([3]) into Eq. ([1]), one obtains the final expression for the 
partition function of a polypeptide experiencing an a-helix<->random coil phase transition. 
This is the expression which we then use to evaluate all thermodynamical characteristics of 
our polypeptide system. 




FIG. 1: Dihedral angles <p and ip used for characterization of the secondary structure of a polypep- 
tide chain. 

B. Molecular dynamics 



Molecular dynamics (MD) is an alternative approach which can be used for the study 
of phase transitions in macromolecular systems. Within the framework of MD, one tries to 
solve the equations of motion for all particles in the system interacting via a given potential. 



Since the technique of MD is well known and described in numerous textbooks 34 
we will only present the basic equations and ideas underlying this method. 
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MD simulations usually imply the numerical solution of the Langevin equation 
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rriiSLi = mjTi 



dU(R) 
<9r ; 



AVi + 7j(t). 



(5) 



Here rrii, Vi and ai are the mass, radius vector, velocity and acceleration of the atom i. 
U(R) is the potential energy of the system. The second term describes the viscous force 
which is proportional to the particle velocity. The proportionality constant $ = 77V7, where 
7 is the damping coefficient. The third term is the noise term that represents the effect of a 
continuous series of collisions of the molecule with the atoms in the medium. To study the 
time-evaluation of the system, the Langevin equations of motion, Eq. (jSJ), are integrated for 
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each particle. 

In this paper, we use the MD approach to study the a-helix<->random coil phase transition 
in alanine polypeptides and compare the results with those obtained using the statistical 
mechanics approach. For the simulations, we use the CHARMM27 force field 39( to describe 
the interactions between atoms. This is a common empirical field for treating polypeptides, 



proteins and lipids 
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MD simulations allow one to study the a-helix<->random coil phase transition of alanine 
polypeptide as this process occurs on the nanosecond time scale. From these simulations, 
one can obtain the important characteristics of the phase transition, such as the transition 
temperature, maximal heat capacity, the temperature range of the transition and the latent 
heat. 

We perform MD simulations of alanine polypeptides consisting of 21, 30, 40, 50 and 100 
amino acids. For this study it is necessary to specify the initial conditions for the system, 



i.e. to define the initial positions of all atoms and set their initia 



velocities. We assume the 



initial structure of the polypeptides as an ideal a-helix [261 . |44j, |45J and assign the particle 
velocities randomly according to the Maxwell distribution at a given temperature. 

The MD simulations of the polypeptides were performed at different temperatures. For 
an alanine polypeptide consisting of 21 amino acids, 71 simulations were performed for the 
temperatures in the region of 300 K° to 1000 K°. For polypeptides consisting of 30, 40, 50 
and 100 amino acids, 31 simulations were performed for each polypeptide in the temperature 
region of 300 K° to 900 K°. The simulations were carried out within a time interval of 100 ns 
and an integration step of 2 fs. The first 25 ns of the simulation were used to equilibrate the 
system, while the next 75 ns were used for obtaining data about the energy and structure 
of the system at a given temperature. 



The set of the parameters used in our simulations can be found in Refs. 34^ 
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All simulations were performed using the NAMD molecular dynamics program 35|], while 



36|. 



visualization of the results was done with VMD 46(. The covalent bonds involving hydrogen 
atoms were considered as rigid. The damping coefficient 7 was set to 5 ps _1 . The simulations 
were performed in the NVT canonical ensemble using a Langevin thermostat with no cutoff 
for electrostatic interactions. 
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III. RESULTS AND DISCUSSION 



In this section we present the results of calculations obtained using our statistical me- 
chanics approach and those from our MD simulations. In subsection IIII Al we discuss the 
accuracy of this force field as applied to alanine polypeptides. In subsection IIII Bl we present 
the PESs for different amino acids in alanine polypeptide calculated versus the twisting de- 
grees of freedom tp and ip (see Fig. [I]). In subsection IIII C[ the statistical mechanics approach 
is used for the description of the a- helix -^random coil phase transition. Here, the results of 
the statistical mechanics approach are compared to those obtained from MD simulations. In 
subsection HHP I the statistical independence of amino acids in the polypeptide is discussed. 



A. Accuracy of the molecular mechanics potential 



The PES of alanine polypeptides was calculated using the CHARMM27 force field [39] 
that has been parameterized for the description of proteins, in particular polypeptides, and 
lipids. Nevertheless, the level of its accuracy when applied to alanine polypeptides cannot 
be taken for granted and has to be investigated. Therefore, we compare the PESs for alanine 
tri- and hexapeptide calculated using the CHARMM27 force field with those calculated using 
ab initio density functional theory (DFT). In the DFT approach, the PES of alanine tri- 
and hexapeptides were calculated as a function of the twisting degrees of freedom, <p and i/j 



(see Fig. [T]), in the central amino acid of the polypeptide 321 ] . All other degrees of freedom 
were frozen. 

To establish the accuracy of the CHARMM27 force field, we have calculated the PESs of 
alanine polypeptides in its /5-sheet conformation. The geometry of alanine tri- and hexapep- 
tide used in the calculations are shown in Fig. [2^ and Fig. [2b respectively. The ab initio 



calculations were performed 32J using B3LYP, Becke's three-parameter gradient-corrected 
exchange functional with the gradient-corrected correlation functional of Lee, Yang and 



Parr 



48j. The wave function of all electrons in the system was expanded using a standard 



basis set B3LYP/6-31G(2d,p). The PESs calculated within the DFT approach have been 



analyzed in Ref. 32]. 



The difference between the PESs calculated with the CHARMM27 force field and with 
the B3LYP functional is shown in Fig. [3] for the alanine tripeptide (left plot) and for the 
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a) 




b) 




FIG. 2: Optimized geometries of alanine polypeptide chains: a) Alanine tripeptide; b) Alanine 
hexapeptide in the /3-sheet conformation. 

alanine hexapeptide (right plot). 

From Fig. [31 we can see that the energy difference between the PESs calculated with the 
CHARMM27 force field and with the B3LYP functional is less than 0.15 eV. To describe 
the relative deviation of the PESs, we introduce the relative error of the two methods as 
follows: 

B3LYp(<P,1p) - EcHARMM27{<P,jj)\d(pdll) _ 
/ \E B 3LYP(P, 1p) + E C HARMM27(P, ^)|d<pd0 

where Eb3lyp(<P,4>) and Echarmm27(<P,'4>) ar e the potential energies calculated within the 
DFT and molecular mechanics methods respectively. Calculating rj for alanine tri- and 
hexapeptide, one obtains: rj 3x Aia = 27.6 % and ?76xA£a = 23.4 % respectively. These values 
show that the molecular mechanics approach is reasonable for a qualitative description of 
the alanine polypeptide. Note however, that the PES obtained for alanine hexapeptide 
within the molecular mechanics method is closer to the PES calculated within the DFT 
approach. This occurs because the PESs Echarmm27('^,4') and Ebslyp^,^) of alanine 
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CHARMM27 B3LYP 

Alanine tripeptide Alanine hexapeptide 
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FIG. 3: Difference between the PESs calculated with the CHARMM27 force field and with the 
B3LYP functional 



321 ] for the alanine tripeptide (left) and the alanine hexapeptide (right). The 
relative energies are given in eV. The equipotential lines are shown for the energies -0.10, -0.05 0, 
0.05 and 0.1 eV. 

hexapeptide were calculated for the structure optimized within the DFT approach, while 
the PESs E C harmm27 and E B3LY p of alanine tripeptide were calculated for the structure 
optimized within the molecular mechanics method and the DFT approach respectively. 

Our analysis shows that the molecular mechanics potential can be used to describe qual- 
itatively the structural and dynamical properties of alanine polypeptides with an error of 
about 20 %. In the present paper, we have calculated the thermodynamical properties of 
alanine polypeptides with the use of MD method and compared the obtained results with 
those attained from the statistical approach. However, ab initio MD calculations of alanine 
polypeptides are hardly possible on the time scales when the a-helix<->random co il p hase 



transition occurs, even for systems consisting of only 4-5 amino acids |30|, |3l|, |32|, |33|, [49] . 
Therefore, we have performed MD simulations for alanine polypeptides using molecular me- 
chanics forcefield. In order to establish the accuracy of the statistical mechanics approach, 
the PES used for the construction of the partition function was also calculated with the 
same method. 
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B. Potential energy surface of alanine polypeptide 

To construct the partition function Eq. ([!]), one needs to calculate the PES of a single 
amino acid in the bounded, e^((p,ip), and unbounded, e^ u '((f,ip), conformations versus the 
twisting degrees of freedom if and ift (see Fig. [p. The potential energies of alanine in different 
conformations determine the and Z u contributions to the partition function, defined in 
Eqs. (T5])-©. 




9 (Degrees) 

FIG. 4: PESs for different amino acids of alanine polypeptide consisting of 21 amino acids calculated 
as the function of twisting dihedral angles ip and ip in: a) second alanine, b) third alanine, c) fourth 
alanine d) fifth alanine and e) tenth alanine. Amino acids are numbered starting from the NH2 
terminal of the polypeptide. Energies are given with respect to the lowest energy minimum of the 
PES in eV. The equipotential lines are shown for the energies 1.8, 1.6, 1.4, 1.2, 1.0, 0.8, 0.6, 0.4 
and 0.2 eV. 

The PES of an alanine depends both on the conformation of the polypeptide and on the 
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amino acid index in the chain. The PES for different amino acids of the 21-residue alanine 
polypeptide calculated as a function of twisting dihedral angles ip and ip are shown in Fig. HI 
These surfaces were calculated with the use of the CHARMM27 forcefield for a polypeptide 
in the a-helix conformation. The PESs a), b), c), d) and e) in Fig. H] correspond to the 
variation of the twisting angles in the second, third, fourth, fifth and tenth amino acids of 
the polypeptide respectively. Amino acids are numbered starting from the NH 2 terminal of 
the polypeptide. We do not present the PES for the amino acids at boundary because the 
angle ip is not defined for it. 

On the PES corresponding to the tenth amino acid in the polypeptide (see Fig. 0fe), one 
can identify a prominent minimum at p> — —81° and ip = —71°. This minimum corresponds 
to the a— helix conformation of the corresponding amino acid, and energetically, the most 
favorable amino acid configuration. In the a— helix conformation the tenth amino acid is 
stabilized by two hydrogen bonds (see Fig. [5]). With the change of the twisting angles 
ip and ip, these hydrogen bonds become broken and the energy of the system increases. 
The tenth alanine can form hydrogen bonds with the neighboring amino acids only in the 
a— helix conformation, because all other amino acids in the polypeptide are in this particular 
conformation. This fact is clearly seen from the corresponding PES Fig. Hfe, where all local 
minima have energies significantly higher than the energy of the global minima (the energy 
difference between the global minimum and a local minimum with the closest energy is 
AE=0.736 eV, which is found at ip = 44° and ip = -124°). 

The PES depends on the amino acid index in the polypeptide. This fact is clearly seen 
from Fig. HI The three boundary amino acids in the polypeptide form a single hydrogen 
bond with their neighbors (see Fig. [5]) and therefore are more weakly bounded than the 
amino acids inside the polypeptide. The change in the twisting angles p and ip in the 
corresponding amino acids leads to the breaking of these bonds, hence increasing the energy 
of the system. However, the boundary amino acids are more flexible then those inside the 
polypeptide chain, and therefore their PES is smoother. 

Fig. H] shows that the PESs calculated for the fourth, fifth and the tenth amino acids are 
very close and have minor deviations from each other. Therefore, the PESs for all amino 
acids in the polypeptide, except the boundary ones can be considered identical. 

Each amino acid inside the polypeptide forms two hydrogen bonds. However since these 
bonds are shared by two amino acids, there is only effectively one hydrogen bond per amino 
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one H-bond *-""^ 

FIG. 5: Alanine polypeptide in the a-helix conformation. Dashed lines show the hydrogen bonds 
in the system. Fig. shows that the second alanine forms only one hydrogen bond, while the fifth 
alanine forms two hydrogen bonds with the neighboring amino acids. 



acid (see Fig. [5j). Therefore, to determine the potential energy surface of a single amino acid 
in the bounded, e^(cp,ijj), and unbounded, e( u '((p,ifj), conformations, we use the potential 
energy surface calculated for the second amino acid of the alanine polypeptide (see Fig. H^i), 
because only this amino acid forms single hydrogen bond with its neighbors (see Fig. [5]). 

The PES of the second amino acid Fig. has a global minima at p = —81° and ip = —66°, 
and corresponds to the bounded conformation of the alanine. Therefore the part of the PES 
in the vicinity of this minima corresponds to the PES of the bounded state of the polypeptide, 
e^ b \(p,tjj). The potential energy of the bounded state is determined by the energy of the 
hydrogen bond, which for an alanine is equal to Ehb =0.142 eV. This value is obtained from 
the difference between the energy of the global minima and the energy of the plateaus at 
V? G (-90°.. - 100°) and ip G (0°..60°) (see Fig. Hi). Thus, the part of the potential energy 
surface which has an energy less then Ehb corresponds to the bounded state of alanine, 
while the part with energy greater then Ehb corresponds to the unbounded state. 

In Fig. [6] we present the potential energy surfaces for alanine in both the bounded (plot 
a) and unbounded (plot b) conformations. Both PESs were calculated from the PES for the 
second amino acid in the polypeptide, which is shown in plot c) of Fig. [61 
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a) oc-helix b) random coil 
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FIG. 6: PESs for alanine in a— helix (plot a) and random coil conformation (plot b). The potential 
energy surface for the second amino acid of the polypeptide is shown in plot c) and is used to 
determine the PESs for alanine in a— helix and random coil conformations. The part of the PES 
shown in plot c, with energy less then Ehb corresponds to the a— helix conformation (bounded 
state) of the alanine, while the part of the potential energy surface with energy greater then Ehb 
corresponds to the random coil conformation (unbounded state). The energies are given in eV. 
The equipotential lines in plot a) are shown for the energies 0.05 and 0.1 and 0.15 eV; in plot b) 
for the energies 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 and 0.9 eV; in plot c) for the energies 1.8, 1.6, 
1.4, 1.2, 1.0, 0.8, 0.6, 0.4 and 0.2 eV. 
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C. a-helix^random coil phase transition in alanine polypeptide 

1. Internal energy of alanine polypeptide 

Knowing the PES for all amino acids in the polypeptide, one can construct the partition 
function of the system using from Eq. (JTJ). Plots a) and b) in Fig. [6] show the dependence 
of e( b '((p,ijj) and €^ u \ip,ip) on the twisting angles ip and ip, while and e- u > define the 
contributions of the bounded and unbounded states of the polypeptide to the partition 
function of the system (see Eqs. (EJ-Q). The expressions for Zi, and Z u are integrated 
numerically and the partition function of the polypeptide is evaluated according to Eq. (CQ). 
The partition function defines all essential thermodynamical characteristics of the system as 
discussed in Ref. [lj]. 

The first order phase transition is characterized by an abrupt change of the internal 
energy of the system with respect to its temperature. In the first order phase transition 
the system either absorbs or releases a fixed amount of energy while the heat capacity as a 
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50]. We study the manifestation 



function of temperature has a pronounced peak |26l . 
of these peculiarities for alanine polypeptide chains of different lengths. 

Fig. [7] shows the dependencies of the internal energy on temperature calculated for alanine 
polypeptides consisting of 21, 30, 40, 50 and 100 amino acids. The thick solid lines corre- 
spond to the results obtained using the statistical approach, while the dots show the results 
of MD simulations. From Fig. UJ it is seen that the internal energy of alanine polypeptide 
rapidly increases in the vicinity of a certain temperature corresponding to the temperature 
of the first order phase transition. The value of the step-like increase of the internal energy 
is usually referred as the the latent heat of the phase transition denoted as Q. The latent 
heat is the energy that the system absorbs at the phase transition. Fig. [7J shows that the 
latent heat increases with the growth of the polypeptide length. This happens because in 
the a-helix state, long polypeptides have more hydrogen bonds than short ones and, for the 
formation of the random coil state, more energy is required. 

The characteristic temperature region of the abrupt change in the internal energy (half- 
wight of the heat capacity peak) characterizes the temperature range of the phase transition. 
We denote this quantity as AT. With the increase of the polypeptide length the dependence 
of the internal energy on temperature becomes steeper and AT decreases. Therefore, the 



16 



phase transition in longer polypeptides is more pronounced. In the following subsection we 
discuss in detail the dependence of AT on the polypeptide length. 
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FIG. 7: Dependencies of the internal energy on temperature calculated for the alanine polypeptide 
chains consisting of 21, 30, 40, 50 and 100 amino acids. Thick solid lines correspond to the results 
obtained within the framework of the statistical model. Dots correspond to the results of MD 
simulations, which are fitted using Eq. (JT]). The fitting functions are shown with thin solid lines. 
The fitting parameters are compiled in Tab. U 

With the molecular dynamics, one can evaluate the dependence of the total energy of the 
system on temperature, which is the sum of the potential, kinetic and vibrational energies. 
Then the heat capacity can be factorized into two terms: one, corresponding to the internal 
dynamics of the polypeptide and the other, to the potential energy of the polypeptide 
conformation. The conformation of the polypeptide influences only the term related to the 
potential energy and the term corresponding to the internal dynamics is assumed to be 
independent of the polypeptides conformation. 

This factorization allows one to distinguish from the total energy the potential energy 



17 



term corresponding to the structural changes of the polypeptide. The formalism of this 
factorization is discussed in detail in Ref. [l|. The energy term corresponding to the internal 
dynamics of the polypeptide neither influence the phase transition of the system, nor does 
it grow linearly with temperature. The term corresponding to the potential energy of the 
polypeptide conformation has a step-like dependence on temperature that occurs at the 
temperature of the phase transition. Since we are interested in the manifestation of the 
phase transition we have subtracted the linear term from the total energy of the system 
and consider only its non-linear part. The slope of the linear term was obtained from the 
dependencies of the total energy on temperature in the range of 300-450 K°, which is far 
beyond the phase transition temperature (see Fig. [7]). Note that the dependence shown in 
Fig. [7] corresponds only to the non-linear potential energy terms. 

The heat capacity of the system is defined as the derivative of the total energy on tem- 
perature. However, as seen from Fig. [7J the MD data is scattered in the vicinity of a certain 
expectation line. Therefore, the direct differentiation of the energy obtained within this ap- 
proach will lead to non-physical fluctuations of the heat capacity. To overcome this difficulty 
we define a fitting function for the total energy of the polypeptide: 



, N AE 
E(T) = Eq H arctan 

7T 



T-Tn 



aT, (7) 

7 J 

where E , AE, T , 7 and a are the fitting parameters. The first and the second terms are 
related to the potential energy of the polypeptide conformation, while the last term describes 
the linear increase of the total energy with temperature. The fitting function E q. (171) was 



used for the description of the total energy of polypeptides in earlier papers [15] , |5l|. The 
results of fitting are shown in Fig. [7J with the thin solid lines. The corresponding fitting 
parameters are compiled in Tab. [H 

Fig. [7J shows that the results obtained using the MD approach are in a reasonable agree- 
ment with the results obtained from the the statistical mechanics formalism. The fitting 
parameter AE corresponds to the latent heat of the phase transition, while the temper- 
ature width of the phase transition is related to the parameter 7. With the increase of 
the polypeptides length, the temperature width of the phase transition decreases (see 7 in 
Tab. [T|, while the latent heat increases (see AE in Tab. H]). These features are correctly 
reproduced in MD and in our statistical mechanics approach. 

Furthermore, MD simulations demonstrate that with an increase of the polypeptide 
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TABLE I: Parameters used in Eq. ([7]) to fit the results of MD simulations. 



11 




LsH/ J 7T 


7 


T. 


a 


21 


11.38±0.24 


1 37±0 10 

X • * J J — 1 — V / a X V / 


79.4±7.6 


670 0±2 

V / | V / a V 1 — I immt a V / 


0471±0 0003 


30 


13.61±0.58 


1.50±0.16 


37.9±7.3 


747.4±3.3 


0.0699±0.0008 


40 


16.80±0.39 


1.991±0.083 


26.6±2.2 


785.7±1.8 


0.0939±0.0005 


50 


19.94±0.79 


2.59±0.21 


29.4±5.5 


786.6±2.9 


0.118±0.0010 


100 


29.95±0.67 


4.00±0.16 


10.5±2.0 


801.1±1.1 


0.2437±0.0009 



length, the temperature of the phase transition shifts towards higher temperatures (see 
Fig- W)- The temperature of the phase transition is described by the fitting parameter T 
in Tab. [B Note also, that the increase of the phase transition temperature is reproduced 
correctly within the framework of the statistical mechanics approach, as seen from Fig. [71 

Nonetheless, the results of MD simulations and the results obtained using the statisti- 
cal mechanics formalism have several discrepancies. As seen from Fig. [7J the latent heat of 
the phase transition for long polypeptides obtained within the framework of the statistical 
approach is higher than that obtained in MD simulations. This happens because within 
the statistical mechanics approach, the potential energy of the polypeptide is underesti- 
mated. Indeed, long polypeptides (consisting of more than 50 amino acids) tend to form 
short-living hydrogen bonds in the random coil conformation. These hydrogen bonds lower 
the potential energy of the polypeptide in the random coil conformation. However, the 
"dynamic" hydrogen-bonds are neglected in the present formalism of the partition function 
construction. 

Additionally, the discrepancies between the two methods arise due to the limited MD 
simulation time and to the small number of different temperatures at which the simulations 
were performed. Indeed, for alanine polypeptide consisting of 100 amino acids 26 simulations 
were performed, while only 3-5 simulations correspond to the phase transition temperature 
region (see Fig. [7J. 
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2. Heat capacity of alanine polypeptide 

The dependence of the heat capacity on temperature for alanine polypeptides of different 
lengths is shown in Fig. [HJ The results obtained using the statistical approach are shown 
with the thick solid line, while the results of MD simulations are shown with the thin solid 
line. Since the classical heat capacity is constant at low temperatures, we subtract out this 
constant value of the for a better analysis of the phase transition in the system. We denote 
the constant contribution to the heat capacity as C300 and calculate it as the heat capacity 
value at 300 K°. The C300 values for alanine polypeptides of different length are compiled 
in the second column of Tab. [Ill 

TABLE II: Parameters, characterizing the heat capacity peak in Fig. [8] calculated using the sta- 
tistical approach. Heat capacity at 300 K, C300, the transition temperature To, the maximal value 
of the heat capacity Co, the temperature range of the phase transition AT and the specific heat Q 
are shown as a function of polypeptide length, n. 



n 


C7 30 o (meV/K) 


To (K) 


Co (eV/K) 


AT (K) 


Q (eV) 


21 


1.951 


740 


0.027 


90 


1.741 


30 


2.725 


780 


0.051 


75 


2.727 


40 


3.584 


805 


0.084 


55 


3.527 


50 


4.443 


815 


0.123 


50 


4.628 


100 


8.740 


835 


0.392 


29 


8.960 



As seen from Fig. [HJ the heat capacity of the system as a function of temperature acquires 
a sharp maximum at a certain temperature corresponding to the temperature of the phase 
transition. The peak in the heat capacity is characterized by the transition temperature To, 
the maximal value of the heat capacity Co, the temperature range of the phase transition 
AT and the latent heat of the phase transition Q. These parameters have been extensively 
discussed in our preceding paper [lj. Within the framework of the two-energy level model 
describing the first order phase transition, it is shown that: 



20 




j i 1 1 " i tt' - u.u r 1 ~ 1 i 

600 700 800 900 750 800 850 900 



Temperature (K°) 



FIG. 8: Dependencies of the heat capacity on temperature calculated for the alanine polypeptides 
consisting of 21, 30, 40, 50 and 100 amino acids. The results obtained using the statistical approach 
are shown with the thick solid line, while the results of MD simulations are shown with the thin 
solid line. Dashed lines show the heat capacity as a function of temperature calculated within the 
framework of the Zimm-Bragg theory 4|. C300 denotes the heat capacity at 300 K°, which are 
compiled in table [III 



AE 

T ~ = const 



C ~ AS 2 ~ n 2 (8) 

Q ~ AE ~ n 
AE 1 



AT 



AS 2 n 

Here AE and AS are the energy and the entropy changes between the a— helix and the ran- 
dom coil states of the polypeptide, while n is the number of amino acids in the polypeptide. 
Fig. [9] shows the dependence of the a-helix^random coil phase transition characteristics on 
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the length of the alanine polypeptide. The maximal heat capacity Cq and the temperature 
range of the phase transition AT are plotted against the squared number of amino acids (n 2 ) 
and the inverse number of amino acids (-) respectively, while the temperature of the phase 
transition To and the latent heat of the phase transition Q are plotted against the number of 
amino acids (n). Squares and triangles represent the phase transition parameters calculated 
using the statistical approach and those obtained from the MD simulations respectively. 

M Statistical model 

A MD approach 




FIG. 9: Phase transition parameters Cq, AT, Tq and Q calculated as a function of polypeptide 
length. Squares and triangles represent the phase transition parameters calculated using the sta- 
tistical approach and those obtained from the MD simulations respectively. 

The results obtained within the framework of the statistical model are in a good agreement 
with the results obtained on the basis of MD simulations. The relative deviation of the 
phase transition parameters calculated in both methods is on the order of 10% for short 
polypeptides and 5% for long polypeptides, as follows from Fig. [9j However, since the MD 
simulations are computationally time demanding it is difficult to simulate phase transition 
in large polypeptides. The difficulties arise due to the large fluctuations which appear in 
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the system at the phase transition temperature and to the large time scale of the phase 
transition process. The relative error of the phase transition temperature obtained on the 
basis of MD approach is in the order of 3 — 5%, while the relative error of the heat capacity 
is about 30% in the vicinity of the phase transition (see Fig. [8]). 

At present, there are no experiments devoted to the study of phase transition of alanine 
polypeptides in vacuo, but such experiments are feasible and are already planned 60]. 
In Ref. ]j| the temperature of the a-helix<->random coil phase transition was calculated. 
Depending on the parameter set, the temperature of the transition ranges from 620 K° to 
650 K° for right-handed a-helix, and from 730 K° to 800 K° for a left-handed a-helix. 

In our previous work 2j on to the theoretical study of phase transitions in polypeptide 
chains, we have introduced the basic ideas of a theoretical method which we have described 
in detail in Ref. [l] and which we currently apply in this work. It was shown that the PES 
calculated as a function of twisting degrees of freedom if and ip determines the partition 
function of the system. To illustrate our method, we used the PES calculated for alanine 

n 

hexapeptide within the framework of the ah initio density functional theory[2|] and obtained 
the phase transition temperature equal to 300 K°. On the other hand, in this paper we 
established that the phase transition temperature of alanine polypeptide in vacuo is 795 K°. 
This is because in Ref. Q the PES was calculated for alanine from the hexapeptide. The 
hydrogen bonds which stabilize the a-helix structure of the hexapeptide are impaired and 
therefore the PES of a single alanine is smoother compared to a long polypeptide where 
every amino acid forms two hydrogen bonds. The smoothing of the potential energy surface 
results in lowering of the energy barriers and the phase transition temperature. 

Nonetheless, smoothing of the PES of an alanine should happen in solution, as the effec- 
tive number of hydrogen bonds in the polypeptide decreases. This fact was demonstrated 
previously [^j], where we compared results of our calculation with available experimental 
data on alanine rich peptides in water solution and observed a good correspondence of the 
phase transition temperature. 

The heat capacity peak is asymmetric. The heat capacity at higher temperatures, beyond 
the heat capacity peak, is not zero and forms a plateau (see Fig. [8). The plateau is formed 
due to the conformations of the amino acids with larger energies 2J. At T=1000 K°), the 
difference in the heat capacity of the polypeptide is 7.6 • 10~ 4 , 1.2 ■ 10~ 3 , 1.6 • 10 -3 , 2.1 ■ 10 -3 
and 4.3 ■ 10~ 3 eV/K° for the Ala2i, Ala 30 , Ala4o, Ala5o and Alaioo peptides respectively. 
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The magnitude of the plateau increases with the growth of the polypeptide length. This 
happens because the number of energy levels with high energies rapidly increases for longer 
polypeptide chains. 



3. Calculation of the Zimm-Bragg parameters 



An alternative theoretical approach for the study of a-helixv->random coil phase transition 
in polypeptides was introduced by Zimm and Bragg 4J. It is based on the construction of the 
partition function of a polypeptide involving two parameters s and a, where s describes the 
contribution of a bounded amino acid relative to that of an unbounded one, and a describes 
the entropy loss caused by the initiation of the a-helix formation. 

The Ztan-Bragg theory Q is se mi e m p iri ea. because it is pa— dependent. The 
theoretical method described in our preceding paper [1[ and which we use in the present 
paper is different as it does not include any parameters and the construction of the partition 
function is based solely on the PES of a polypeptide. Therefore, the construction of our 
partition function is free of any parameters, and this is what makes it different from the 
models suggested previously. Assuming that the polypeptide has a single helical region, the 
partition function derived within the Zimm-Bragg theory, reads as: 



n— 3 



Q = T + a J2(n-k-2)s k 



(9) 



k=l 



where n + 1 is the number amino acids in the polypeptide, s and a are the parameters of 
the Zimm-Bragg theory. The partition function, which we use in the present paper Eq. (pQ) 
can be rewritten in a similar form: 



Z 



(n-l)-3 

l + ps(T) 3 (n-k-3)s(Tf 

k=l 



(10) 



Here n is the number of amino acids in the polypeptide and the functions s(T) and £,{T) 
are defined as: 
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s(T) 



ten 



/TV f'TX 
— 7T J — 7T 



kT 



d^d-0 



/7T P7T 
— 7T J —IT 



exp 



e( u )((^,i/>) 

w 



kT 



d^d^ 



(11) 



(12) 



where (cp, ip) and tp) are the potential energies of a single amino acid in the bounded 

and unbounded conformations respectively calculated versus its twisting degrees of freedom 
tp and ip. By comparing Eqs. Q and ffTUl) . one can evaluate the Zimm-Bragg parameters as: 



a(T) = (3(T)s(T) 3 , (13) 

where j3(T) is defined in Eq. (j4]). 

The dependence of the Zimm-Bragg parameters s and a on temperature is shown in 
Fig. [TUh and Fig. [TUb respectively. The function — i?Tln(s) grows linearly with an increase 
in temperature, as seen in Fig. [TUh . The zero of this function corresponds to the temperature 
of the phase transition in an infinitely long polypeptide. In our calculation it is 860 K° (see 
black line in Fig. [TUh). Parameter a is shown in the logarithmic scale and has a maximum 
at T = 560 K°. Note, that this maximum does not correspond to the temperature of the 
phase transition. 



The parameters of the Zimm-Bragg theory were considered in earlier papers 16, 
In Fig. [TUh we present the dependence of parameter s on temperature calculated in 
(see squares, triangles and stars in Fig. [TUb) using a matrix approach described in Ref. 
The energies of different polypeptide conformations were calculated using the force field 
described in Ref. 53j. Squares, triangles and stars correspond to three different force field 




parameter sets used in Ref. [ijj] , which are denoted as sets A, B and C. Fig. [TUh shows that 
the results of our calculations are closer to the results obtained using the parameter set C. 
This figure also illustrates that the Zimm-Bragg parameter s depends on the parameter set 
used. Therefore, the discrepancies between our calculation and the calculation performed in 
Ref. jl9| arise due to the utilization of different force fields. 

The Zimm-Bragg parameter a was also calculated in Ref. 19|. However, it was not 



systematically studied for the broad range of temperatures, and therefore we do not plot it 
in Fig. [TUb. In Ref. 19[] the parameter a was calculated only for the temperature of the 
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a-helix<-i>random coil phase transition ranging from 620 K° to 800 K°. In Ref. 19j, it was 
also demonstrated that parameter a is very sensitive to the force field parameters, being in 
the range 10 -9 — 10~ 3 6 . In our calculation a = 10~ 3 4 at 860 K°. The dependence of the 
parameter a on the force field parameters was extensively discussed in Ref. 19|] , where it was 
demonstrated that this parameter does not have a strong influence on the thermodynamical 
characteristics of phase transition. 
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FIG. 10: Dependence of the parameters of the Zimm-Bragg theory [4| s (plot a) and a (plot b) on 
temperature. Parameter s describes the contribution to the partition function of a bounded amino 
acid relative to that of an unbounded one. The parameter a describes the entropy loss caused by 
the initiation of the o-helix formation. Parameter s was also calculated in Ref. [r| using three 
different force fields, shown with stars, triangles and squares in plot a. 



If the parameters s and a are known, it is possible to construct the partition function of 
the polypeptide in the form suggested by Zimm and Bragg [l], and on its basis calculate all 
essential thermodynamic characteristics of the system. The dependence of the heat capacity 
calculated within the framework of the Zimm-Bragg theory is shown in Fig. [8] by dashed 
lines for polypeptides of different length. 
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From Fig. [HJ it is seen that results obtained on the basis of the Zimm-Bragg theory are 
in a perfect agreement with the results of our statistical approach. The values of the phase 
transition temperature and of the maximal heat capacity in both cases are close. The 
comparison shows that the heat capacity obtained within the framework of the Zimm-Bragg 
model at temperatures beyond the phase transition window is slightly lower than the heat 
capacity calculated within the framework of our statistical model. 

An important difference of the Zimm-Bragg theory from our theory arises due to the 
accounting for the states of the polypeptide with more than one a— helix fragment. These 
states are often referred to as multihelical states of the polypeptide. However, their statistical 
weight in the partition function is suppressed. The suppression arises because of entropy 
loss in the boundary amino acids of a helical fragment. The boundary amino acids have 
weaker hydrogen bonds than amino acids in the central part of the a-helix. At the same 
time the entropy of such amino acids is smaller than the entropy of an amino acids in the 
coil state. These two factors lead to the decrease of the statistical weight of the multihelical 
states. 

The contribution of the multihelical states to the partition function leads to the broaden- 
ing of the heat capacity peak while the maximal heat capacity decreases. The multihelical 
states become important in longer polypeptide chains that consist of more than 100 amino 
acids. As seen from Fig. [HJ the maximal heat capacity obtained within the framework of the 
Zimm-Bragg model for Alaioo polypeptide is 10% lower than that obtained using our sug- 
gested statistical approach. For alanine polypeptide consisting of less than 50 amino acids 
the multihelical states of the polypeptide can be neglected as seen from the comparison per- 
formed in Fig. [HJ Omission of the multihelical states significantly simplifies the construction 
and evaluation of the partition function. 



4- Helicity of alanine polypeptides 



Helicity 
imentally 



is an important characteristic of the polypeptide which can be measured exper- 
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24|. It describes the fraction of amino acids in the polypeptide that 



are in the a-helix conformation. With the increase of temperature the fraction of amino 
acids being in the a— helix conformation decreases due to the a-helix^random coil phase 
transition. In our approach, the helicity of a polypeptide is defined as follows: 
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^n- 4/ 



where n is the number of amino acids in the polypeptide, Z&, Z u are the contributions to 
the partition function from amino acids in the bounded and unbounded states defined in 
Eqs. (T5|) and ([3j) respectively. The dependence of helicity on temperature obtained using the 
statistical approach for alanine polypeptides of different length are shown in Fig. [TTJ 
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FIG. 11: Dependency of the helicity on temperature obtained using the statistical approach for 
alanine polypeptide chains consisting of 21, 30, 40, 50 and 100 amino acids. The helicity for alanine 
polypeptide consisting of 21 amino acids obtained within a framework of MD approach is shown 
in the inset. 

On the basis of MD simulations, it possible to evaluate the dependence of helicity on 
temperature. Helicity can be defined as the ratio of amino acids being in the a-helix con- 
formation to the total number of amino acids in the polypeptide, averaged over the MD 
trajectory. The amino acid is considered to be in the conformation of an a-helix if the an- 
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gles describing its twisting are within the range of </? G [—72°; — 6°] and ip G [0°; —82°]. This 
region was chosen from the analysis of angles <p and if) distribution at 300 K°. The helicity 
for alanine polypeptide consisting of 21 amino acids obtained within the framework of MD 
approach is shown in the inset to Fig. [TTJ From this plot it is seen that at T f» 300 K°, which 
is far beyond the temperature of the phase transition, the helicity of the Ala2i polypeptide 
is 0.82. The fact that at low temperatures the helicity of the polypeptide obtained within 
the MD approach is smaller than unity arises due to the difficulty of defining the a-helix 
state of an amino acid. Thus, the helicity obtained within the MD approach rolls off at 
lower temperatures compared to the helicity of the polypeptide of the same length obtained 
using the statistical mechanics approach. 

The kink in the helicity curve corresponds to the temperature of the phase transition of 
the system. As seen from Fig. [TTJ with an increase of the polypeptide length, the helicity 
curve is becomes steeper as the phase transition is getting sharper. In the limiting case of 
an infinitely long polypeptide chain, the helicity should behave like a step function. This is 
yet another feature of a first-order phase transition. 



An important question concerns the statistical independence of amino acids in the 
polypeptide at different temperatures. In the present section we analyze how a particular 
conformation of one amino acids influences the PES of other amino acids in the polypeptide. 
In Fig. [12] we present the deviations of angles ip and if> from the twisting angles cpio and ipio 
in the 10 — th amino acid of alanine polypeptide. These results were obtained on the basis 
of MD simulations of the Ala2i polypeptide at 300 K° and at 1000 K°. The deviation of 
angles ip and if) is defined as follows: 



3=1 

where i is the amino acid index in the polypeptide and M is the number of MD simulation 
steps. Note, that the plots shown in Fig. [12] do not depend on the reference amino acid (we 



D. Correlation of different amino acids in the polypeptide 




(14) 
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used the middle amino acid in the polypeptide). 
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FIG. 12: The root mean square deviation of angles ip and tp calculated with the use of Eq. (TTj 
for alanine polypeptide consisting of 21 amino acids. The calculations were done in respect to the 
tenth amino acid of the polypeptide for 300 K (top plot) and for 1000 K (bottom plot). 

The top plot in Fig. [Ti2l was obtained at 300 K°. At this temperature, all amino acids 
in the polypeptide are in the a— helix conformation, and the deviation of angles <p and ip is 
less than 16° for all amino acids except the boundary ones, where the relative deviation of 
the angles tp and ip is 28° and 34° respectively. This happens because, while the boundary 
amino acids are loosely bounded, the central amino acids in the polypeptide are close to 
the minima that corresponds to an a— helix conformation. In the a— helix state, all central 
amino acids are stabilized by two hydrogen bonds, while the boundary amino acids form 
only one hydrogen bond. 

At 1000 K° the polypeptide is, to large extent, found in the random coil phase and 
therefore becomes more flexible. In the random coil phase, the stabilizing hydrogen bonds 
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are broken, and the deviation of angles ip and ip significantly increases. This fact is clearly 
seen in the bottom plot of Fig. [121 However at 1000 K, the deviation of angles <p and ip in 
the central and in the boundary amino acids is almost the same, confirming the assumption 
that in the random coil phase, short alanine polypeptides do not build hydrogen bonds. 

Another important fact which is worth mentioning is that in the random coil phase (and 
in the central part of the a— helix), the deviation of angles <p and ip does not depend on 
the distance between amino acids in the polypeptide chain. For instance, the deviation 
between angles in the 10 — th and in the 11 — th amino acid is almost the same as the 
deviation between angles in the 10 — th and in the 17 — th amino acid. This fact allows one 
to conclude that in a certain phase of the polypeptide (a-helix or random coil), amino acids 
can be treated as statistically independent. 



IV. CONCLUSION 



In the present paper we presented results of calculations obtained with the statistical 
method described in our preceding paper l|. We have also performed a detail analysis 
of the a-helix<->random coil transition in alanine polypeptides of different lengths. We 
have calculated the potential energy surfaces of polypeptides with respect to their twisting 
degrees of freedom and constructed a parameter-free partition function of the polypeptide 
using our statistical formalism jl|] . From this partition function, we derived and analyzed the 
temperature dependence of the heat capacity, latent heat and helicity of alanine polypeptides 
consisting of 21, 30, 40, 50 and 100 amino acids. Alternatively, we have obtained the 
same thermodynamical characteristics from the use of molecular dynamics simulations and 
compared them with the results of our statistical mechanics approach. The comparison 
proved the validity of our method and established its accuracy. 

It was demonstrated that the heat capacity of alanine polypeptides has a peak at a certain 
temperature. The parameters of this peak (i.e. the maximal value of the heat capacity, the 
temperature of the peak, the width at half maximum, the area of the peak) were analyzed 
as a function of polypeptide length. Based on the predictions of the two energy-level model, 
it was demonstrated that the a-helix ^random coil transition in alanine polypeptide is a 
first order phase transition. 

We have established a correspondence of our method with the results of the semiempir- 
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ical approach suggested by Zimm and Bragg [4]. For this purpose we have determined the 
key parameters of the Zimm-Bragg semiempirical statistical theory The calculated param- 
eters of the Zimm-Bragg theory were compared with the results of earlier calculations from 
Ref. [19]. 

The final part of this paper deals with the statistical independence of amino acids in the 
polypeptide at different temperatures. It was shown that a particular conformation of one 
amino acids influences the PES of other amino acids in the polypeptide. We demonstrated 
that in a certain phase, amino acids can be treated as statistically independent. 

In this paper, we demonstrated that the new statistical approach is applicable for the 
description of a-helix<-i>random coil phase transition in alanine polypeptides. However, this 
method is general and can be used to study similar processes in other complex molecular 
systems. For example, it would be interesting to apply the suggested formalism to the study 
of /?-sheet<->random coil phase transition and to the study of non-homogeneous polypeptides 
(i.e. consisting of different amino acids). The suggested method can also be applied to the 
description of protein folding — an important question left open for further consideration. 

In this work we have investigated a-helix^->random coil phase transition of alanine 
polypeptides in vacuo. So far there has been done no experimental work on a-helix<->-random 

;hat such experiments are possible 
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57] and the ESI mass 



coil transition in gas phase. Nevertheless, it is important 
and can be performed using of the techniques MALDI 
spectroscopy 58|, |59J . We hope that our theoretical analysis of the a- helix -^random coil in 
alanine polypeptides in vacuo will stimulate experimentalists to verify our predictions. 
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